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AUTOMATED METHOD TO REMOVE FALSE ECHOES DUE TO SOLAR 

INTERFERENCE 



FIELD OF THE INVENTION 

[0001] The present invention relates to weather radar processing and display. In 
particular, the present invention is directed to a method of removing false echoes from weather 
radar data caused by solar interference. 

BACKGROUND OF THE INVENTION 

[0002] As is known in the art, NEXRAD radar sites perform elevation scans at various 
angles to collect weather-related data. Fig. 1 illustrates the various elevation angles that a 
NEXRAD radar site may scan and the range for each scan angle. NEXRAD radar sites may 
experience interference caused by the radar detecting electromagnetic radiation from another 
source, most often another radar or the Sun. The Sun particularly poses problems near sunrise 
and sunset at elevation scan angles approximately equal to the angle of the Sun. The interference 
generally corrupts one or more radials of radar data, and particularly affects most or all of a 
radial, usually out to the maximum range. Such interference from the Sun is known in the art as 
a "Death Ray," which will be used herein to mean false echoes that occur near sunrise and sunset 
which are caused by solar radiation. Examples of the false echoes are shown in Figs. 2 and 3, as 
"Death Rays" 10, 12, 14, 16 and 18. 

[0003] A problem is that NEXRAD radar data provided by the National Weather 
Service does not contain any logic to detect and/or correct for interference from the Sim. The 
problem is so large that over 300 of the NEXRAD 0.5 degree Base Reflectivity products are 
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affected by the Death Ray effect each day. The number increases if all NEXRAD elevation 
scans are included. This problem results in degraded radar products, as is evidenced by Figs. 2 
and 3, where the inclusion of the Death Rays 10, 12, 14, 16 and 18 makes the product more 
difficult to interpret. 

[0004] Another problem is that radials are wider as they go out in range (e.g., the 
radials look similar to a baseball bat, with the narrow end towards the radar and the fat end 
pointing at the source of contamination). Thus, it is likely that real echoes will be embedded 
within the false echoes. When interference is encountered, the reflectivity values increase with 
range, with a maximum being about 35 dBZ. The interference appears for one or two volume 
scans for a single elevation and up to 30 minutes on a Composite Reflectivity product. Other 
reflectivity derived products, such as echo tops, layer reflectivity and composite reflectivity are 
also affected. Still further, radial velocity and spectrum width products usually indicate range- 
aliased data, and due to radial velocity corruption, the mesoscale process can be negatively 
affected. 

[0005] Yet another problem is that many of the NEXRAD site clocks are not 
synchronized to a standard time, thus making the calculation of the Sun's position less accurate. 
Thus, without a certain time reference, automatic identification and removal of the Death Ray is 
difficult as the Sun's position relative to the particular radar scan cannot be accurately 
determined. 

[0006] Given the right tools and time, a trained Meteorologist can manually identify 
and remove the false echoes. This is the approached used by providers such as WSI of Billerica, 
MA for their NOWrad product which is available every 15 minutes. However, there is a 
significant drawback to this method as it delays the release of the product for up to 7 minutes. 
This is a significant amount of time, especially when there is fast changing weather. As a result, 
for more timely products, such as WSI's 5-Minute NOWrad PLUS and Premium NOWrad, 
manual removal is not practical. 

[0007] Thus, there is a need for a system of automatic removal of Death Ray effects 
that can be quickly and automatically performed. In addition, there is a need for a system of 
automatic removal of Death Ray effects that is very accurate and does not remove real echoes. 
The present invention provides for such a system. 

SUMMARY OF THE INVENTION 

[0008] The present invention is directed to systems and methods for automated 
detection and removal of Death Rays due to solar interference. In accordance with the present 
invention, the removal process is applied in real time to NEXRAD BREF1 products or other 
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similar radar products. The removal process does not impact the timeliness of critical products, 
and conservatively removes solar interference. In other words, the process may leave in some 
false echoes rather than to remove real echoes. Still further, the present invention overcomes 
problems with inaccuracies with NEXRAD clock time synchronization. 

[0009] In accordance with an aspect of the present invention, there is provided a 
method of removing solar interference from radar data. The method includes extracting a site 
latitude, a site longitude and scan elevation time from the radar data; determining the position of 
the Sun for the extracted latitude and longitude and scan time; determining a radial that has been 
contaminated with solar interference; and removing the solar interference. 

[0010] In accordance with a feature of the invention, the removal process is only 
performed if the Sun's position is within a threshold elevation angle of an elevation scan angle. 
The threshold elevation angle may be between 0.5 and 1.3 degrees. 

[0011] In accordance with another feature, the method accounts for inaccuracies in a 
clock at a radar site from which the radar data was received. This may include identifying the 
radial closest to a Sun azimuth and its neighbors and sorting a predetermined number of radials 
closest to the Sun position based on a number of non-zero echoes in each radial. 

[0012] In accordance with yet another feature, determining a radial contaminated with 
solar interference may include determining if the number of non-zero echoes in a highest priority 
radial exceeds a second predetermined threshold; determining if a difference of the number of 
non-zero echoes between the highest and lowest radial exceed a percentage threshold; and 
determining if a three highest priority radials are azimuthally consecutive, and if so, returning a 
middle radial, otherwise returning a radial with the highest number of non-zero echoes. 

[0013] In accordance with another feature, removing the solar interference further may 
include for each range gate of a radial N identified as being contaminated by solar interference, 
determining if an echo in a current range gate exceeds an intensity threshold; and for echoes in 
the current range gate that do not exceed the intensity threshold, examining the current range 
gates for radials N-2 and N+2 determine if an echo for the current range gate for those radials is 
zero. The intensity threshold may be variable over a length and increase as the range increases 
for a radial. If the range gates for the radials N-2 and N+2 are not zero, then the value of the 
current range gate for radial N is not changed and the next range gate for the radial N is 
examined. However, if the range gates for the radials N-2 and N+2 are zero, then the value of 
the current range gate for radial N is changed to zero and the next range gate for the radial N is 
examined. 
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[0014] In addition, for each range gate of the N+l and N-l radials, the process may 
determine if an echo in a current range gate in the N+l or N-l radial does not exceed the 
intensity threshold; and if an echo in the current range gate in the N+l or N-l radial does not 
exceed the intensity threshold, examine if the current range gates for both radials N-l and N+l 
are zero; and set a value of the current range gate for the radial N to 0 if the current range gates 
for both radials N-l and N+l are zero. 

[0015] In accordance with another aspect of the invention, there is provided a method 
of determining if radar data contains solar interference. The method includes extracting a site 
latitude, a scan longitude and scan elevation time from the radar data; determining the position of 
the Sun for the extracted latitude and longitude and scan time; determining if an elevation angle 
of the Sun is within a threshold elevation angle of a scan elevation angle; and examining radar 
data for radials closest to the Sun's position if the Sun is within the threshold elevation angle. 

[0016] In accordance with a feature of the invention, the threshold elevation angle is 
between 0.5 and 1.3 degrees. 

[0017] In accordance with another feature, examining radar data for radials closest to 
the Sun's position further may include determining if the number of non-zero echoes in a highest 
priority radial exceeds a second predetermined threshold; determining if a difference of the 
number of non-zero echoes between the highest and lowest radial exceed a percentage threshold; 
and determining if a three highest priority radials are azimuthally consecutive, and if so, 
returning a middle radial, otherwise returning a radial with the highest number of non-zero 
echoes. 

[0018] In accordance with yet another feature, the method may further include 
removing the solar interference from radials containing the solar interference in accordance with 
the first aspect of the invention. 

[0019] Additional features and advantages of the invention will be made apparent from 
the following detailed description of illustrative embodiments that proceeds with reference to the 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0020] The foregoing summary, as well as the following detailed description of 
preferred embodiments, is better understood when read in conjunction with the appended 
drawings. For the purpose of illustrating the invention, there is shown in the drawings exemplary 
constructions of the invention where like elements have like reference numerals; however, the 
invention is not limited to the specific methods and instrumentalities disclosed. In the drawings: 
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[0021] Fig. 1 is a graph of elevation angles vs. range for NEXRAD radar sites; 

[0022] Figs 2 and 3 are exemplary radar products showing Death Rays; 

[0023] Figs. 4 and 5 are the exemplary radar products of Figs. 2 and 3, respectively, 
after the solar interference removal process of the present invention has removed the Death Rays; 

[0024] Figs. 6-8 are flow charts illustrating the processes performed to remove Death 
Rays from radar data; and 

[0025] Fig. 9 illustrates an exemplary radar installation, radials and range gates. 

DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS 

[0026] The solar interference removal process will now be described with reference to 
Figs. 6-9. The process may be performed in any combination of computer hardware and/or 
software, as would be recognized by one of ordinary skill in the art. Referring now to Fig. 6, 
there is a high level flow chart of the process to remove solar interference from NEXRAD Level 
III Data. At step 100 the process begins by receiving NEXRAD Level III BREF I Data. At step 
102, the NEXRAD site latitude and longitude and elevation scan time is extracted. At step 104, 
the Sun azimuth and elevation for the input latitude and longitude and time is determined. The 
determination of the Sun's azimuth is well known by those of ordinary skill in the art, and can be 
calculated based on start time of elevation scan and site latitude/longitude. For example, the 
Sun's position can be determined using the Solar Position Calculator at the National Oceanic and 
Atmospheric Administration Website. After the Sun's position is determined for the input 
latitude/longitude and time, it is determined if the Sun elevation is within a threshold elevation 
angle (ELEV). In accordance with the present invention, the ELEV threshold is between 0.5 and 
1.3 degrees. In other words, the removal process is not run if the Sun elevation is more then 0.5 
to 1.3 degrees of the scan elevation angle. Preferably, the value of ELEV is 0.5, however, it may 
need to be larger to compensate for situations where the NEXRAD clock is not correct. Thus, at 
step 106, if the Sun elevation is outside 0.5 to 1.3 of the scan elevation angle, then the process 
stops at step 108 as it is unlikely that solar interference has occurred. If the Sun's position is 
within the ELEV threshold, then the interference removal procedure (part 1) is run at step 1 10 to 
identify the radial having the greatest interference. 

[0027] The procedure run at step 110 will now be explained in greater detail with 
reference to Fig. 7. The process of step 1 10 is preferably run because of the inaccuracies of 
NEXRAD clocks. At steps 120 and 122, from the NEXRAD Level III data, the radial closest to 
the Sun azimuth is identified and its neighbors. At step 124, the 7 radials closest to the computed 
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Sun position are sorted based on the number of non-zero echoes in each radial. It is noted that a 
different number of radials may be sorted at step 124. 

[0028] At step 126, it is determined If the number of echoes in the highest priority 
radial exceeds a SIZE threshold (e.g., 20). If not, then the routine exits at step 128 and returns 
"no interference detected." If the number of echoes exceeds the SIZE threshold, then it is 
determined if the difference of the number of echoes between the highest and lowest radial 
exceeds a DIFFERENCE percentage threshold (e.g., 10%). If not, then the routine exits at step 
132 and returns "no interference detected." If the difference of the number of echoes between 
the highest and lowest radial exceeds the DIFFERENCE threshold at step 103, then it is 
determined at step 134 if the three highest priority radials are azimuthally consecutive, and if so, 
the azimuth of the middle radial is returned at step 138. In accordance with the present 
invention, for a three radial contamination, the solar removal process identifies the middle radial. 
Otherwise, if the three highest priority radials are not consecutive, then at step 136 the azimuth 
of radial with the highest number of non-zero echoes is returned. 

[0029] The processing returns to step 112, where it is determined if the process of Fig. 
7 returned an interference pattern (i.e., either the azimuth of the highest priority radial or the 
middle radial of the three azimuthally consecutive highest priority radials). If so, then the 
interference removal procedure (part 2) is run at step 116, otherwise, the process ends at step 
114. 

[0030] The interference removal procedure of step 1 16 is shown in greater detail with 
reference to Figs. 8 and 9, where from the NEXRAD Level III data (step 140) a two phase 
procedure is applied to the radials having interference. At step 142, for each range gate of the 
radial(s) identified as having interference, it is determined at step 144 if an echo in the current 
range gate exceeds an INTENSITY threshold. The INTENSITY threshold is variable over the 
length and increases as the range increases. The INTENSITY threshold is set at 25 dBZ when 
the range is less then 230 kM. Outside of 230 kM, the INTENSITY threshold is gradually 
increased to a maximum value (e.g. 40 dBZ) at the furthest range. The INTENSITY threshold 
provides a rough estimate of the presence of solar radiation as the Sun provides a constant signal 
which results in a greater computed dBZ value the further you go out in range. It is noted that 
the intensity of a Death Ray rarely exceeds 35 dBZ. 

[0031] If at step 144 the echo in the current range gate exceeds the INTESITY 
threshold, the value is not changed and the next range gate is examined. If the echo in the 
current range gate does not exceed the INTESITY threshold, then at step 146, the current range 
gates for radials N-2 and N+2 (see, Fig. 9, where N is the current radial being examined) are 
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examined to see if the echo for the current range gate for those radials is zero. It is assumed that 
if solar interference is present for radial N, that it will likely only affect radial N and radials N +/- 
1, but not radials +/-2. Thus, if the range gates for the radials N-2 and N+2 are not zero, then the 
value of the current range gate for radial N is not changed at step 152 and the next range gate for 
the radial is examined at step 142. Otherwise, at step 150 the value of the current range gate for 
radial N is set to 0 at step 150 and the next range gate for the radial is examined at step 142. 

[0032] Once all of the range gates are examined at steps 142-152, then the phase 2 is 
performed beginning at step 154 where for each range gate of the N+l and N-l radial(s), it is 
determined at step 156 if an echo in the current range gate exceeds an INTENSITY threshold. 
The INTENSITY threshold is the same as noted above with regard to step 144. If at step 156 the 
echo in the current range gate exceeds the INTESITY threshold, the value is not changed and the 
next range gate is examined at step 154. If the echo in the current range gate does not exceed the 
INTESITY threshold, then at step 158, the current range gates for radials N-l and N+l are 
examined to see if the echo for the current range gate for those radials is zero. Thus, if the range 
gates for the radials N-l and N+l are not zero, then the value of the current range gate for radial 
N is not changed at step 164 and the next range gate for the radial is examined at step 142. 
Otherwise, at step 162 the value of the current range gate for the radial is set to 0 at step 162 and 
the next range gate for the radial is examined at step 154. Once all range gates are examined at 
steps 154-164, the process exits at step 166 to return to step 118, where the main routine exits. 

[0033] In accordance with the present invention, an examination of a radial and its 
neighbors is performed such that only the radial data that actually was affected by solar radiation 
is removed. This way, the Death Ray is removed while true echoes remain. The operation of the 
present invention can been seen by examining Figs. 2 vs. 4 and Figs. 3 vs. 5, where the removal 
of death rays 10, 12, 14, 16 and 18 is evident, while the real echoes are retained. 

[0034] While the present invention has been described in connection with the preferred 
embodiments of the various Figs., it is to be understood that other similar embodiments may be 
used or modifications and additions may be made to the described embodiment for performing 
the same function of the present invention without deviating therefrom. For example the process 
could be applied to other NEXRAD Level III elevation scans and products (i.e., velocity), 
NEXRAD Level II data, or other radars (i.e. TDWR, customer radars, etc.). 

[0035] Furthermore, it should be emphasized that a variety of computer platforms, 
including handheld device operating systems and other application specific operating systems are 
contemplated. Still further, the present invention may be implemented in or across a plurality of 
processing chips or devices, and storage may similarly be effected across a plurality of devices. 
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Therefore, the present invention should not be limited to any single embodiment, but rather 
should be construed in breadth and scope in accordance with the appended claims. 
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